function l = l_f(pm,x,e,p,w,t,r)

l= (pm.rho.*(1-t).*p./wf(pm,x,e,w,r).*x.*e).^(1./(1-pm.rho));